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Abstract Solutions to Einstein's field equations describing rotating fluid bod- 
ies in equilibrium permit parametric (i.e. quasi-stationary) transitions to the 
extreme Kerr solution (outside the horizon) . This has been shown analytically 
for discs of dust and numerically for ring solutions with various equations of 
state. From the exterior point of view, this transition can be interpreted as a 
(quasi) black hole limit. All gravitational multipole moments assume precisely 
the values of an extremal Kerr black hole in the limit. In the present paper, 
the way in which the black hole limit is approached is investigated in more 
detail by means of a parametric Taylor series expansion of the exact solution 
describing a rigidly rotating disc of dust. Combined with numerical calcula- 
tions for ring solutions our results indicate an interesting universal behaviour 
\Q , of the multipole moments near the black hole limit. 
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1 Introduction 

A fascinating property of some equilibrium configurations of rotating fluid 
bodies is the existence of a "black hole limit" when a well-defined parameter 
relation is reached. This was first demonstrated by Bardeen and Wagoner [1] 
with their approximate solution to the problem of a uniformly rotating disc 
of dust in general relativity. The exact solution to the disc problem found by 
Neugebauer and Meinel [2,3] confirmed this result, see [4,5] for details. In the 
limit, the spacetime separates: From the "exterior point of view" , the extreme 
Kerr metric outside the event horizon is formed. From the "interior point of 



A. Kleinwachter • H. Labranche ■ R. Meinel 
Theoretisch-Physikalisches Institut, University of Jena, 
Max-Wien-Platz 1, 07743 Jena, Germany 
E-mail: meinel@tpi.uni-jena.de 



view" , a regular, non-asymptotically flat spacetime - containing the fluid body 
- with the extreme Kerr "near-horizon geometry" (also called "throat geome- 
try" [6]) at spatial infinity results. Similar phenomena were observed for lim- 
iting solutions to the static Einstein- Yang-Mills-Higgs and Einstein-Maxwell 
equations [7—11]. Strictly speaking, there is not yet a horizon in the limit, 
and hence the denotation "quasi-black hole" used by Lemos et al. [9] is also 
appropriate. However, as already discussed by Bardeen [12], the slightest dy- 
namical perturbation will lead to a genuine black hole. Therefore, we continue 
to use the term "black hole limit" . Such a limit was also found numerically for 
uniformly rotating fluid rings with various equations of state [13-15]. 

It was shown in [16] that the black hole limit of a (uniformly) rotating fluid 
body in equilibrium always leads to an extremal Kerr black hole 1 . A necessary 
and sufficient condition for the limit is the parameter relation [18] 

M - 2flJ ->• 0, (1) 

where M denotes the total (gravitational) mass, J the angular momentum and 
Q the angular velocity 2 . Note that Lemos and Zaslavskii [19] have shown in 
a more general context that quasi-black holes are always extremal if infinite 
surface stresses are excluded. 

In [15] evidence for a universal behaviour of the Geroch-Hansen gravita- 
tional multipole moments [20,21] near the black hole limit of rotating fluid 
bodies was reported. In the present paper, these results are confirmed and 
extended. In particular, for the disc of dust, a parametric Taylor series ex- 
pansion of the exact solution at the black hole limit is derived. The multipole 
moments can be calculated according to the algorithm introduced by Fodor 
et al. [22] using the Ernst formulation of the stationary and axisymmetric 
vacuum Einstein equations [23,24]. 



2 The Ernst potential of the rigidly rotating disc of dust 

2.1 The metric and the Ernst potential 

The asymptotically flat spacetime of a rigidly rotating disc of dust surrounded 
by a vacuum can be described using the metric 

ds 2 = e- 2t/ [e 2 W + dC 2 ) + pW] - c 2U (ad V + dt) 2 , (2) 

where the functions e 2fe , e 2U and a depend only on the canonical Weyl coor- 
dinates p and C, and not on the time coordinate t and the azimuthal angle if, 
corresponding to stationarity and axial symmetry The equatorial "plane" is 
given by £ = and the axis of rotation by p = 0. 



See [17] for a discussion of the extension of the Kerr uniqueness proof to the case with 
a degenerate horizon. 

2 Throughout the paper we use units in which the speed of light as well as Newton's 
gravitational constant arc equal to 1. 
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Fig. 1: The thick line is the infinitesimally thin disc of dust. The disc has a 
radius po an d rotates about the £-axis. 



Einstein's vacuum field equations, in the case of stationarity and axial 
symmetry, are equivalent to the complex Ernst equation [23,24] 

(R/)V 2 / = V/ • V/ , (3) 

where f(p, £) is a complex function called the Ernst potential, and V 2 and V 
are respectively the Laplace and the gradient operators in a three dimensional 
Euclidean space, as if p, Q and ip were cylindrical coordinates. The metric 
potentials can then be calculated from: 

e 2U = K/ , 
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where b{p, Q is the imaginary part of the Ernst potential b := 9/ (so / 
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-\b). For asymptotically flat solutions with the gravitational mass M and 



angular momentum J, the behaviour of the Ernst potential at infinity is 
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where r := -\/p 2 + C 2 an d tan0 := p/C,. 



2.2 The rigidly rotating disc of dust 

The thin disc of dust is represented in Fig. 1. It is located at ( = 0, p < po and 
rotates with the constant angular velocity J? with respect to infinity about the 
C-axis. It is useful to introduce a potential V by 



-,2V 



= c 2U [(l + f2a) 2 -f2 2 p 2 e~ 



For any uniformly rotating fluid body in equilibrium 3 the function V is con- 
stant (V = Vb) along the surface of the body (the disc "consists" of a surface 
only) . This constant Vb is related to the relative redshift Zq of photons (with- 
out ^-component of orbital angular momentum) emitted from the surface and 
received at infinity via 

c- Va -1 = Z Q . (10) 

The surface mass density is not constant within the disc: it must vary in a 
unique way such that each particle of dust can follow a geodesic of constant p, 
all with the same Q : by the sole effect of gravitation. The physical parameters 
that we just mentioned can be combined into a "relativistic parameter" called 
p: 

p:=2Q 2 plc- 2Vo . (11) 

The disc has a set of solutions that can be specified by fixing two parame- 
ters. But by using normalized and dimcnsionless coordinates, the disc can be 
characterized by a single parameter. Thus, from its Ernst potential we can 
extract dimcnsionless combinations of physical parameters as functions of p 
alone. The gravitational mass M, the angular momentum J, the imaginary 
part 60 of the Ernst potential at the centre of the disc (p — 0, ( — + ) and 
the previously introduced parameters are given by 

9V h'cn 2 (I,h') , N 

C 2V„ = ±^_i ^ ( 12a ) 

sn(/,//)dn(/,/i') 
°o = 7 . ( 12b ) 



1 / h' 2 

Q Q = Qp = -^\- — cn(/, ti) , (12c) 

M = 2QM = -b - i? ci , (12d) 

J=A[2 2 J = -b -2J? ci , (12c) 

where sn(u, k), cn(u,k) and dn(u, k) are the Jacobian elliptic functions, and 
h, h! , I and c\ are the following functions of ji: 



h' = \M-^==), (12f) 





In Vl + x 2 + x) 

9{x) = . ; 1 , 12g 

^Jp-x Vl + x 2 

ci= -^{2^/1 + p 2 E(am(J, ti),h') -{p+^l + p 2 )I + /1} (12h) 

1 f g(x)x n , N 

and I n = - - , v ; Ht . (12i) 

* Jo V M - x 



3 While for the disc of dust the metric (2) can be used globally, it can only be introduced 
in the exterior (vacuum) region in case of genuine fluid bodies. 



The Ernst potential for a uniformly rotating disc of dust is given in [25,26] for 
the whole spacetime in terms of hyperelliptic theta functions and quadratures. 
Since the multipole moments of stationary and axisymmetric configurations - 
which contain all the information about the exterior gravitational field - can 
be obtained from the Ernst potential on the symmetry axis, it is sufficient to 
restrict ourselves to this axis, where the potential becomes simpler, see [27]. 
The representation of the potential that we will use is considerably different 
from the one given in [27], but all the necessary relations are provided in detail 
in [28,5]. 

To make the Ernst potential of the disc a function of one physical param- 
eter, in our case p, we begin by introducing the normalized and dimensionless 
coordinates x := pj po and y := C/Po- The potential on the axis is obtained by 
evaluating it in the limit x — > 0; it then becomes a function of two variables: 
f(p;x = 0,y) — > f(p;y). By restricting ourselves to the positive part of the 
axis (y > 0), the Ernst potential can be written as follows: 
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with the following functions defined as real functions of p and y: 
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Y = sign(y - t)F(v, ti) , 
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and i?2 being the Jacobian theta function defined by 



-d 2 (x;B) = ^^ ex P ' 



-(2»+l) 



B + (2n+l)x 



The disc solution is physically relevant from the Newtonian limit given by 
fi — > to the ultra-relativistic limit, i.e. the black hole limit, given by the 
smallest positive value of \i for which cn(7, h') — 0. This "upper limit" of /z is 
called /io, i.e. 

< H < Mo = 4.6296618434743420427 . . . (14) 



2.3 The black hole limit of the disc 

In the black hole limit, the source shrinks to the origin of the Weyl coordinate 
system, i.e. for the disc we have po — ^ 0. Therefore, it is appropriate to intro- 
duce two new pairs of dimcnsionlcss Weyl coordinates in addition to x and y, 
one normalized with twice the angular velocity fi and a second one normalized 
with the mass M (note that 2fiM — >• 1 in the black hole limit, hence the two 
normalizations become identical in the limit): 



p:=2J?p, C := 2f?C and 
9 ' M ' M ' 



(15a) 
(15b) 



We will also use this "tilde" and "hat" notation to make some physical param- 
eters dimcnsionlcss. The logic of this notation is that a "tilde" introduces pow- 
ers of 2J7 such that the parameter becomes dimcnsionless (mass M = 2J?M, 
angular momentum J = 4J? 2 J, etc.), and the "hat" notation does the same 
with powers of M (mass M = 1, angular momentum J = J/M 2 etc.). 

These normalized coordinates can be introduced in the Ernst potential of 
Eq. (13) by rewriting y with the help of 2J?M = M(n) and Qp = f2 (/j,) in 
one of the following two ways: 



A) 

C_ 

A) 



2QM C 



M(p) 



2Qp M 2Qq{h) 

2(2( C 



C, 



2f2p a 



2^o (ji) 



(16a) 
(16b) 
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The black hole limit of the disc is reached for p — > hq, or equivalently for 
c Va — > 0. As discussed in the Introduction, there is a separation of spacetime 
in the limit. We are interested here in the "outer world" given by the extreme 
Kerr metric outside the horizon. Note that the horizon of the extreme Kerr 
metric is situated at p = ( = 0. On the axis, the outer world corresponds to 
finite values of (, or £, while the "inner world" is described by finite values of 
y. If we evaluate all the functions from Eq. (13a) in the black hole limit for 
C > 0, we find: 

N = l, Q±= 1±C- (17) 

This gives a potential which is identical to the Ernst potential of the extremal 
Kerr black hole on the (upper part of the) axis: 

/(WC) = |— ^- (18) 

C + i-i 

This potential can uniquely be extended off the axis and reads 

r — 1 — i cos 9 
f = ?+l-icos0 ' (19) 



where r := yp 2 + C 2 an d tan# := p/Q. 



2.4 Leading order behaviour of the solution close to the black hole limit 

Taking the Ernst potential of any asymptotically flat solution at p — 0, one 
can expand it at £ — > +oo in the following way: 

1 - /(C) _ ^ rn n 



E 



i + /(C) f^ n+1 ' 

The coefficients m n are the same as those used in [22] to compute the multipolc 
moments defined by Geroch [20] and Hansen [21]. The first two coefficients are 
too = M and m\ — i J . Let us rewrite the coefficients in a normalized and 
dimensionless form: 

with k„ '■= M 



k n * \M 



,71-1 



m„ 



Note that the k n correspond exactly to the coefficients m n of a Kerr solution 
with mass M and angular momentum J, meaning that to* = 1 for black holes. 
It turns out that the following relations hold in the black hole limit: 

= , (20a) 



dp 

d2 * 
TO„ 



dp 2 



= (20b) 

A«=Mo 



or, equivalently, 

dm! 



d 2 



de v ° 

* 

rn„ 



d(e^) 



Vb^ 



= , (21a) 

= . (21b) 

Vo^ — oo 



Actually, numerical evidence suggests strongly that (21) holds for a large class 
of solutions (not only the disc) with parametric transition to a black hole 4 . 
Indeed, Eq. (21a) was verified in [15] (in a different but equivalent form) for 
several sequences of numerical solutions describing rotating fluid rings reaching 
the black hole limit. We have checked that the second property (21b) holds for 
these sequences as well. Note that the relations (21) also hold for the normal- 
ized Geroch-Hanscn multipolc moments (m* replaced by the Geroch-Hansen 
moments normalized with the k n ). This follows directly from the structure of 
the expressions given in [22]. 

It can easily be concluded that the Ernst potential of any fluid configuration 
having a black hole limit (e Vo — > 0) and fulfilling the conditions (21) can be 
decomposed into an Ernst potential /Kerr of the Kerr solution with the same 
mass and angular momentum and a residual potential R of order e 3Vo : 

/(M, J; C) - /Kcrr(M, J; C) + R(M, J; C) , (22a) 

f Km (M,J;Q= ^(C-Mj-iJ ( 22b ) 

J y ' ,s; M(C + M)-iJ ' V ' 

R(M,J;C)= O(e 3V0 ) . (22c) 

Note that /xorr is not restricted to J < M 2 as it is for black holes. Indeed, 
rings and discs have J > M 2 , thus /Kerr is given here by the axis potential of a 
"hyperextreme" Kerr solution. In the above equations one may take M as one 
parameter and consider J as a function of M and of the second parameter e^ , 
i.e. J = J(M, e v °). For the disc solution, J/M 2 is a function of e v ° alone. At the 
black hole limit, c v ° = 0, J = M 2 and R vanishes. Comparing Eqs (22) with 
Eq. (9) one can see that R must also decay rapidly in the far field, R = 0(£~ 3 ), 
corresponding to the fact that /Kerr already contributes the correct first two 
multipolc moments. Eqs (22) can be extended off the axis with /Korr as the 
Ernst potential of the general Kerr solution. While / and /Kerr are then both 
solutions to the Ernst equation, R is of course not - due to the non-linearity. 



3 Taylor series at the black hole limit 

Our aim is to expand the Ernst potential /(/i; C) or /(m; C) °f the disc solution 
at the black hole limit (/i — > /xq) in terms of powers of (/J, — Hq). Such an 



4 It was shown in [18] that e v o — ¥ is necessary and sufficient for reaching a black hole 
limit. An equivalent condition is the relation (1). 
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Fig. 2: Four exact functions of jjl compared to their Taylor approximations 
with orders n = 5 and 10. On the right side of the plot, the functions are from 
the top to the bottom: 2QM, Qpo, e 2V ° and &o- 



expansion can then easily be transformed to an expansion in terms of the 
more universal parameter e v ° . 

To find an explicit form of the Taylor series, we make use of the computer 
algebra systems Maple and Mathematica and we divide the work in two steps. 
First, we write down the series of all functions that depend on /j, alone. Then, 
we do the same for the remaining functions which depend on /i and (, or £. 



3.1 Series of functions of /z 

The set of functions that depend on /z alone is listed in Eqs (12). For each of 
these functions, the resulting derivatives of fi in the limit \x — >• \xq give pure 
numbers. 

The direct computation of a series for I with the form given previously is 
impracticable because zeros come out for some denominators. To avoid this 
pathology, we rewrite the function with the substitution x = [i sin <f>. This also 
has the effect of removing the /i-dcpendence of the upper limit of integration: 



%) 



<A 



9J X ) 
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g(/j,sm 2 (j))sm(f)d(j). (23) 
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-1.1979704X10" 1 
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6.1997318xl0~ 3 


2.8702661x10" 
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8.2373333 xl0~ 3 


-2.2207483xl0~ 2 
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-2.1917290xl0~ 3 


-3.9472326x10" 


-3 


4.3533289 xl0~ 4 


1.1246071 xl0~ 3 
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4.5766410xl0~ 4 


3.5824068x10" 


-4 


-1.9230828xl0~ 4 


1.7087311 xlO~ 5 
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-7.5851829X10- 5 


-2.0388433x10" 


-5 


3.8239383X10- 5 


-1. 4784817X10" 5 
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1.1139165X10- 8 


-9.6924305x10" 


-7 


-6. 5882003x10-° 


3.0205335xl0~ 6 
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-1.5243204xl0~ 6 


6.0947891x10" 


-7 


1.0947125xl0~ 6 


-5.0848591xl0~ 7 


8 


1.9941471X10" 7 


-1.4423676x10" 


-7 


-1.8180382X10" 7 


8.0708022xl0~ 8 


9 


-2.5296376xl0~ 8 


2.7737508x10" 


-8 


3.0627438 xl0~ 8 


-1.2616141X10" 8 


10 


3. 1381072 xl0~ 9 


-4.9011518x10" 


-9 


-5.2656087X10- 9 


1.9805535X10- 9 



Table 1: First coefficients aj of the expansions defined in Eq.(25) for the func- 
tions e 2Vo , b , f}p Q and 2QM. 



The series reads 



I(fi) = 1.5752 + 0.24046 (/i - Mo) - 0.017245 (p 
+ 0.0017270 (m - mo) 3 + C[(n - Mo) 4 ] , 



Mo) 



(24) 



where decimal numbers are truncated after 5 significant digits (we omit the 
ellipsis after each number from here on). Similarly, the derivatives of all other 
functions of m that we need for the Ernst potential can be calculated recur- 
sively. Thus we can already present Taylor series of a few relevant parame- 
ter functions such as the Ernst potential at the centre of the disc, given by 
f(p = 0, C = + ) = /o = c 2V " + i&o, as well as the dimensionless products ilp 
and 1QM. These functions need to be expanded with high numerical precision 
since they enter into the remaining series expansions. 

Let us now introduce the following notation for the n-th order Taylor ap- 
proximation A n of a function of \x near the black hole limit: 



A n (bo(fi)) 



3=0 



Oj(M-Mo)^ 



(25) 



where 6o(m) is used as an example. For the short list of functions that we have 
introduced, we show their first expansion coefficients aj in Table 1. The same 
functions are also plotted in Fig. 2 with their respective Taylor approximations 
of order n = 5 and n = 10. This figure allows us to assess the quality of the 
Taylor series. Approximations with n — 5 and n = 10 are indistinguishable 
from the exact function on the plot for \i > 2, while polynomials of higher 
orders (n = 10 versus n — 5 on the plot) improve the approximation near 
the Newtonian limit (m — > 0) reasonably well. In Fig. 3, one can see that the 
Taylor scries seem to converge to their respective exact functions from the 
black hole limit (m —> Mo) all the way down to the Newtonian limit (fj, —¥ 0). 
The convergence is readily seen in the figure with e 2Vo , while for the three other 
functions, one must take into account the logarithmic scale to appreciate it. 
Since the function e 2V ° is related to the redshift parameter as given in Eq. (10), 
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we can calculate how wrong the redshift becomes from the approximations. 
In the Newtonian limit, where the Taylor approximations have the greatest 
deviations, the redshifts, approximated by Zq = [*A n (e 2Vi, )]~ 1//2 — 1, become 
Z = 0.17691, 9.0567 x 10~ 3 , 4.2302 x 10~ 4 for n = 5, 10, 15 respectively, 
while of course Zq — > is the correct result. 



3.2 Series of functions of /i and £ 

The next step to obtain a series of the disc's axis potential is to expand the 
remaining terms which depend on both /i and £. To this end, we introduce 
the normalized coordinate £ = 2QQ in every function given in Eqs (13). By 
computing series at ^ = no for L(/j,,Q, S(p,(), H(p,(), Z(n, C), A(nX) an d 
T(/z,£), we can then determine the series of N(fi, £), which can be combined 
with the series of Q±(n,C) to obtain the Ernst potential of the disc on the 
axis. 

For the function L((i, £), it i s helpful to carry out the substitution x = 
/xsin 2 <f>, as in Eq. (23), in order to calculate it numerically: 

4^oC 7 2 (C 2 +4^ o 2 )M 2 sin 2 + 4l2 2 . 2 . 

— g(/ism 0) sin dip, 



7r V / (C 2 +4f2 2 )V + 16/# / C 2 + 4r? 2 cos 



where i?o is a function of /i already expanded before. Once the integrand is 
expanded in a series in (/i — /io), the integral becomes easier to perform on 
each individual term of the expansion. This series contains polynomials of l/£ 
with odd exponents and reads 

T . ~, 0.24697 , . 0.015102, , 2 

MM,C) = ~ — (M-Mo) ~ (A* -Mo) 

0.0063269 0.0021138 \ . , 3 mu , 4l 

~ + ~ (^-^of + OKM-Mo) 4 ] , (26) 

C C J / 

where numbers are again truncated after five significant digits. The next func- 
tions to expand, S(fi,C) and -P(/x, £), do not require anything other than a 
direct computation of the series from the computer. The series for S reads 

S(ij,,()= 0.053082- 0.011080 {»- Ho) (27) 

0.0022735- °- 00 ! 2 1646 )(,-,q) 2 

-4.5822 x 10- 4 + °- 00 ! 1646 ) ( M - mo) 3 + 0[(m - Mo) 4 ] • 

C 2 / 

The function H{jjl 1 £) is not difficult to expand by calling P(fx, C) where needed. 
Note that the sign from the factor "sign(y — t)" has to be chosen to be positive 
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1 



(a) \e 2V o - An(e 2Vo )\ 



(b) \b -A n (bo)\ 




(c) \Qp -An(Opo)\ 



(d) \2f2M - A n (2f2M)\ 



Fig. 3: These figures show the deviations of the n-th order Taylor approxima- 
tions An from their respective exact functions. The Taylor series are truncated 
to order n — 1 ... 15 for the functions e 2V ° and bo, and to n = ... 15 for f2p$ 
and 2f2M. The lines are ordered from the smallest n to the largest when 
one follows the abcisssa from right to left. The parameter /i is shown on the 
abscissae and the ordinates show on a logarithmic scale the absolute value 
of the deviations. The plunges to zero for a few curves only mean that the 
corresponding polynomials A n intersect the exact function at that point. 
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since y — > oo for finite £ at the black hole limit. Of course, our series will diverge 
in a certain region near the disc. This is related to the previously discussed 
"separation of spacetime" as /j, — > mo- The result for H is 

0.0013618\ 

T 



H{n, C) = - 0.055908 + ( 0.012264 + - — - I (m - Mo) 



r>^r.,i~ 6.6713 x 10- 4 0.0064947V , 2 
-0.0026412 ^ + ^ (fj, - mo) 

C C 2 / 



5.5869 x 10 



-4 



2.1410 x 10~ 4 0.0023506 



c e 

2.3723 x 10~ 4S 



c 3 



(M-Mo)' + 0[(M-Mo) 4 ]. (28) 



The remaining functions Z(/i, £), A(/j,, £) and T(/i, £) all depend on the elliptic 
integrals given by X(fi, (,) and Y(m, C)- Since we restrict ourselves from now 
on to the above mentioned "plus sign" , the latter functions become simply 
X = F,(v,h') and Y — F(v,h'). A good strategy to produce series of these 
functions is given by expanding the right hand sides of 



v = am(y, h ) = arcsin 
sin(u) = sn(Y, h!) = 
cos(u) = cn(y, h') = 



VI -hP 



hi 
VI -hP 

h' 
y/h(P-h) 



h> 

^Jl-h' 2 sm 2 (v) = dn(Y, h') = VhP . 

Then, the Taylor series of F(v, h') and E(t>, h') can be computed as scries 
containing derivatives of the Jacobian elliptic functions. Each time that such 
derivatives need to be evaluated, the answer can easily be picked up in the 
four series obtained from above. By properly combining the series of F(v, h') 1 
E(v,h'), K(h), K(h') and E(/i'), we obtain: 

Z(n,Q= 0.0028257+ (-0.0011830+ L9350 ^ < 10 j (^ _ /iQ ) 

3.6763 x 10- 4 - L762( L X 10 " 6 - 3 - 3 ° 02 ^ 10 ' 4 ) (m - Mo) 2 

-1.0047 x 10- 4 + 9 - 5M8 * 10 ' 7 + L859 l X 10 ' 4 

c c 2 

;US: ' (,/Uri \M-Mo) 3 + 0[(M-Mo) 4 ], (29) 



C 3 
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., ~, 0.066452. . 0.0085785. , 2 

A(p,C)= ~ (M-Mo) + ~ (m~Mo) 2 (30) 

7.1026 x 10- 4 0.0012715\ , ., mu . 4 . 
~ + ~ (M - Mo) 3 + 0[( P - mo) 4 ] , 

C C 3 / 

„, ~, 0.20935, . 0.0042719,. . 2 

T{H,Q= ~ — (M-Mo) ~ (M-Mo) 

0.0038072 . 0.0040061\ , x3 4 



C C 3 / 



( M - Mo )' + 0[( M - Mo ) 4 ]. (31) 



Note that we obtained the series of T by the calculation of a series of dT/d/z 
using several relations for derivatives of the theta function $2 and subsequent 
integration. At this stage, the remaining computations are straightfoward. The 
series of iV(/x, C) is obtained by combining together the series from Eqs (24) 
and (26-31): 

0.25127, N , ( 0.012990 . 0.031568\ , n2 



N(ji,0= 1-— -^(/z-/z )+ — -~ + — -~ (M-Mo) 

c v c c / 

0.0056452 0.0032641 2.1819 x 10~ 4 \ . 

(M-Mo) 3 



c c 2 c 3 

+ O[( M - M0 ) 4 ]. (32) 

And the expansions of J?o(m) and 60 (/z) are needed for Q±(/z, £), which gives: 

Q±(M.O - l±C-0.028702(zi- Ato ) 2 +0.0039472( M - Ato ) 3 +C[(At-Mo) 4 ] • (33) 

3.3 Series of the Ernst potential of the disc 

Both real series (32) and (33) can now be combined into the complex Ernst 
potential given by Eq. (13a). The series takes the following form: 

~ C-l-i 0-25127C 

/(m;C)= ~ ;--- — — M-Mo 

C + l-i (C+l-i) 2 



+ 



0.063137i 0.012845 + 0.012990 i 0.044414 



(C + l-i) 3 (C + l-i) 2 C + l-i, 
0.015864 -0.0065282 + 0.0064556 i 



(m - Mo) 2 



(C + l-i) 4 (C + l-i) 

0.0034840- 0.022452 i 



(C + l-i 



1 2 



-0.0066618+ 0.0033960 i 
== 0.0022492 

C + l-i 

\4 



(M -^Mo) 3 

C 



0[(/z - mo) 4 ] • (34) 
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(a) e 2U at C = 10 



(b) b at C = 10 




(c) e 2u at ( = 2 



(d) fe at C = 2 




o 1 



(e) e 2u at C = 0.5 



(f) 6 at C = 0.5 



Fig. 4: Real and imaginary parts of the Ernst potential and their respective 
Taylor approximation as function of jjl for ( =10, 2 and 0.5. The solid lines 
are from the exact Ernst potential, while the dotted, dash-dotted and dashed 
lines are Taylor polynomials with orders n — 3, 6, 9 respectively. 
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We only wrote down the beginning of our results and with only five signifi- 
cant digits, since the space needed for further orders inflates rapidly. But the 
method that we used to compute the Ernst potential allows us to generate the 
series beyond ten orders with more then ten significant digits in a reasonable 
amount of time; e.g. a personal computer with a 2.2 GHz CPU takes around 
2 minutes to compute all series up to ten orders and ten significant digits. In 
Fig. 4, the series of the Ernst potential is shown as Taylor polynomials for n = 
3, 6 and 9 at different positions on the axis, and it is compared to the exact 
potential. One can see that for large ^-distances from the disc (in Fig. 4, £ = 
10 or 2), the series seems to converge for any value of /i. Closer to the disc 
(C = 0.5 in Fig. 4), the series does not converge any more for values of /i too 
much smaller than no- One can also see that for very relativistic discs (/i > 2), 
the series with orders like n = 3 or 6 give excellent approximations as long as 
the series is evaluated for sufficiently large £. 

The Taylor series using the other normalized and dimensionless coordinate, 
( = (/M, reads as follows: 

t( ?\ C-l-i 0-0258361 2 

/(M;C)= rTT -- ifTT -^(M-,o) (35) 

-0.0062737+ 0.0032261 if,.. A 3 -.,, , 4l 



C(C + l-i) 2 



-( M - Mo )^ + 0[( M -Mor 



Another useful representation is related to Eqs (22) - taking into consid- 
eration that e Va = 0[(fi — /io)]: 

/Gu;C) = /Kerr(M;C) + R(j*;Q (36) 

with 



/K( 


srrWC) = 


C + i- 


ij(u) 




R(jm;0 = 


-0.0062737 

C(C + l-i) 2 


Note that 








J( M ): 


= 1 + 0.012918 (n- 


Mo) 2 



J( M ) 



-bp(/j) - 2J? (/i)ci(/i) 
[bo(fJ-) + n (fJ.) ci(ix)] 



3 i /-or/.. .. \4i 



{H - no) 6 + 0[(fi - /io 



0.0016130 (n - no) 3 + 0[(p - no 



,4l 



The three representations of the series of the disc's Ernst potential given 
above show that the solution near the limit differs from the extreme Kerr 
solution with the mass Mbh = 1/2/2 by a term of order (/i — /io) [see (34)], 
from the extreme Kerr solution with the same mass M by a term of order 
(/i — no) 2 [ scc (35)] and from the Kerr solution with the same mass M and 
angular momentum J by a term of order (// — no) 3 [ see (36)]. 
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4 Discussion 

The expansion of the axis potential /(/z;C) of the disc solution as given in 
(35) has a remarkably simple structure: the contributions at each order of 
(p — fio) are given by rational functions of £. However, the most interesting 
aspect of our results is the fact that the solution near the limit is approximated 
very well by the (hyperextreme) Kerr solution with the same M and J as the 
disc - up to a residual term of order C(e 3y °), see (22) and (36). A precise 
formulation of this statement is that all multipole moments differ from those 
of the Kerr solution with the same M and J only by terms of order 0(e 3V °). 
The fact that the Kerr black hole, which is the unique stationary black hole 
surrounded by an asymtotically flat vacuum region, is completely characterized 
by the two parameters M and J is called the "no-hair theorem" . Accordingly, 
all multipole moments of the Kerr spacetime are unique functions of M and 
J. When approaching the black hole limit of the disc, the deviations of the 
multipole moments from those of the corresponding Kerr spacetime (same M 
and J) decay oc e 3V ° . Thus one can say, in a sense, that the "hair" vanishes in 
this way. As already stressed in subsection 2.4 this seems to hold for all fluid 
bodies in equilibrium that permit a black hole limit. We also mention that, 
according to our experience, the exterior spacetime of a uniformly rotating 
fluid body is never described by the Kerr metric - except for the black hole 
limit (e Va — > 0). As long as e v ° 7^ 0, all "higher" multipole moments (beyond 
M and J) always seem to be greater than those of the corresponding Kerr 
spacetime [29]. 
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